ApJ submitted, Sept. 14, 2012 

Preprint typeset using I^T^X style cmulateapj v. 5/2/11 



A FAST FLARE AND DIRECT REDSHIFT CONSTRAINT IN FAR-ULTRAVIOLET SPECTRA OF THE 

BLAZAR S5 0716+714 

Charles W. Danforth, Krzysztof Nalewajko, Kevin France, & Brian A. Keeney 

CASA, Department of Astrophysical & Planetary Sciences, University of Colorado, 389-UCB, Boulder, CO 80309; 

danforth@casa.colorado.edu 
ApJ submitted, Sept. 14, 2012 

ABSTRACT 

The BLLacertae object S5 0716+714 is one of the most studied blazars on the sky due to its 
active variability and brightness in many bands, including VHE gamma rays. We present here two 
serendipitous results from recent far-ultraviolet spectroscopic observations by the Cosmic Origins 
Spectrograph onboard the Hubble Space Telescope. First, the blazar increased in flux rapidly by ~ 40% 
(—0.45 magh -1 ) followed by a slower decline (+0.36 mag h" 1 ) to previous far-UV flux levels during 
the course of our 7.3 hour HST observations. We model this flare using asymetric flare templates and 
constrain the physical size, and energetics of the emitting region. Furthermore, the spectral index of 
the object softens considerably during the course of the flare from a„ — 1.0 to a v « —1.4. Second, 
we constrain the source redshift directly using the > 30 intervening absorption systems. A system at 
z = 0.2315 is detected in Lya, Ly/3, O VI, N V, and C III and defines the lower bound on the source 
redshift. No absorbers are seen in the remaining spectral coverage (0.2315 < ZL ya < 0.47) and we 
set a statistical upper bound of z < 0.304 (90% confidence) on the blazar. This is the first direct 
redshift limit for this object and is consistent with literature estimates of z = 0.31 ± 0.08 based on 
the detection of a host galaxy. 

Subject headings: BL Lacertae objects: individual: S5 0716+714, radiation mechanisms: nonthermal, 
galaxies: jets, galaxies: active, intergalactic medium, quasars: absorption lines, 
ultraviolet: general 



1. INTRODUCTION 

Blazars are a class of AGN in which non-thermal con- 
tinuum emission is relativistically beamed into our line 
of sight, rendering them extremely bright across many 
wavebands at cosmological distances. S5 0716+714 is 
one of the most extensively studied blazars in the sky. 
It belongs to the blazar sub-class of BL Lacertae Objects 
(BL Lacs) which are known for their smooth non-thermal 
continua without clear emission lines. This lack of iden- 
tifying spectral features makes the source redshifts, and 
hence the physical characteristics of these astrophysically 
interesting objects, very challenging to determine unam- 
biguously. 

Spectroscopi cally, S5 0716+714 was first studied in the 
optical band bv lBiermann et al~l ()1981h who found a com- 
pletely featureless spectrum. No emission or absorption 
lines were ever identified in this source, despite several 
attempts dStickel et al.lll993t iVermeulen fc Tav"ior]|1995l: 
iRector fc Stockd 120011: iFinke et al.ll2008t) . iNilsson et al.l 
(2008) detected an excess in the point spread func- 
tion of S5 0716+714 which is consistent with a host 
galaxy at z = 0.31 ± 0.08. This range is consistent 
with the iStickel et al.l (|1993t ) identification of a pair of 
galaxies near the sight line at z ~ 0.26, suggesting 
that the host is part of a galaxy group. Detection of 
S5 0716+714 in Very- High-Energy (VHE) gamma-rays 
by the MAGIC telescope led t o different upper limit s 
on the source redshi ft: z < 0.5 (lAnderhub et al.l l2009h , 
and z < 0.21 ± 0.09 (jPrandini et al.ll2010D . though these 
depend on the poorly-characterized cosmic infrared back- 
ground. 

Blazars are also known for variability. Even amongst 



blazars, S5 0716+714 is well known for its strong and 
persis tent optical variabil ity with a duty cycle of ~ 
90% (jOstorero et all 120061) . Intraday variability (IDV) 
is routinely measured in this source with the reported 
maximum variabi lity rates of ~ 0.1 — 0.16 magni- 
tudes per hour (iWaener et all [1991 IWu et al.l 120051: 
i Villata et all l200tt IMontagni et al.l 120061: lOstorero et al.l 
l2006t INesci et all I2002D to as high as 0.38 magnitudes 
per hour ()Chandra et al.l 120111 ) . An extensive statisti- 
cal analysis bv iMontagni et ail (j2006l ) showed that vari- 
ability rates of order ~ O.lmaghr -1 usually last for no 
more than two hours, and their duty cycle is only a 
few percent. Asymmetric flares, with the increase rate 
on average fa s ter th an th e decline ra t e, wer e found by 
IWaener et all (1 19961) and INesci et al.l (I2002D, se e how - 
ever. iGhisellini et al.l (|1997l ) and IMontagni et al.l (|2006t) . 
Quasi-periodic oscillations were claimed over short time 
intervals with various periods: ~ 1 day and ~ 7 days 
([Quirrenbach et alJl!991t iWagner et aTlll996[ ). ~ 4 days 
(IHeidt fc Wagnedll996D . and even as fast as ~ 15 min- 
utes (|Rani et al.ll2010D . 

BLLac objects are extremely useful as probes of in- 
tervening material. High signal-to- noise (S/N) optical 
and ultraviolet spectral observations are possible if the 
blazar is observed during a bright state. Their smooth, 
power-law continua make it possible to distinguish weak, 
broad intergalactic absorption features from weak fea- 
tures in the background source continuum. It was for 
this reason that S5 0716+714 was observed by the Cos- 
mic Origin Spectrograph on th e Hubble Space Telescope 
fHST/COS ^Green et al.ll20T2h . 

The far-ultraviolet (FUV) spectrum of S5 0716+714 
provides an excellent dataset relevant to IGM cosmology 
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Fig. 1.— Coadded COS spectrum of S5 0716+714. Flux has 
been scaled to the first exposure (however, see discussion in Sec- 
tion 3.1) and smoothed by 35 pixels (approximately five resolution 
elements). Shaded data shows coadded photon noise. Sharp emis- 
sion features are geocoronal airglow. The spectral ranges of the 
two medium-resolution, far-UV gratings on COS are shown. 

which will be presented as part of a much larger study 
of the low-redshift IGM (Danforth et al. 2013, in prep). 
We present here two serendipitous discoveries from the 
S5 0716+714 observations. First, our five-orbit HST ob- 
servations happened to observe a fast, asymmetric flare 
during which the blazar brightened by > 35% and then 
faded to approximately the flux level at the start of our 
observations. We discuss the flare and changes to the 
spectral energy distribution (SED) in Section 3. Sec- 
ond, in Section 4, we use the observed intervening ab- 
sorption systems (the Lya forest) along the sight line to 
place the first direct constraints on the source redshift of 
S5 0716+714. We discuss and summarize our results in 
Section 5. 

2. OBSERVATIONS AND TIMING ANALYSIS 

S5 0716+714 was observed with the Cosmic Origins 
Spectrograph (COS) on December 27, 2011, as part of 
HST program 12025 (PI: Green). Two exposures were 
made during each of five HST orbits over the course 
of 7.3 hours; the first five exposures with the G130M 
(1135 < A < 1450 A, 6.00 ksec) grating, the second 
five with the G160M (1400 < A < 1795 A, 8.25 ksec) 
gratings. The ten exposures were obtained from the 
Mikulski Archive for Space Telescopes (MAST) and re- 
processed with a recent version of CALCOS (2.17.3a). 
The calibrated, one-dimensional spectra were next coad- 
ded wi th the custom LpL p rocedures described in de- 
tail by iDanforth et al.l (|2010H . We present the coadded 
COS spectrum of S5 0716+714 in Figure [T] (however, see 
the cautions regarding coaddition of a strongly varying 
source in Section 3.1.) 

A Voigt profile fit to the Galactic Lya profile in the 
G130M data gives a column density of log N(HI) — 
20.372 ± 0.003. This corresponds to a color excess of 
E( B — V) = . 045 an d we deredden all observed fluxes 
via iFitznatrickl (flQOWl . 

Many AGN show some degree of variability between 
one epoch and another and our coaddition software au- 
tomatically scales exposures from different epochs to cor- 
rect this (usually minor) effect. However, during coad- 
dition of the S5 0716+714 data, it was noticed that the 

Note that the absorption-derived color excess measured here 
is ~ 50% larger than the extinction value de rived from the low- 
resolution Schlcgcl, Finkbincr, & Davis (1998) dust maps. 



mean flux varied by as much as 35% during the course 
of 7.3 hours. COS is a photon-counting instrument; the 
arrival time of each photon is recorded as well as the 
two-dimensional position on the detector. It is there- 
fore possible to ext ract very high reso lution (~ 30 ms) 
time-domain data (France et al.l 120101 ) ■ We extract an 
[xi,yi,ti] photon list from each exposure i and coadd 
these to create a master photon list. The to- 

tal number of counts in a [Ax, Ay] box is integrated 
over a timestep At. We use a timestep of At = 120 s 
for S5 0716+714 to balance time-resolution with rea- 
sonable photometric precision. Continuum count rates 
are extracted over ~ 30 (37) A integration regions for 
the G130M (G160M) grating at four nominal far-UV 
wavelengths (1190 = [1178.6 - 1208.5 A] and 1390 = 
[1377.8 - 1407.7 A] in G130M; 1460 = [1447.9 - 1484.6 
A] and 1520 = [1507.3 - 1544.0 A] in G160M). Narrow 
interstellar and intergalactic absorption lines are present 
in the data (see Figure [1]) , however these features only 
remove a small number of photons relative to the bright 
continuum emission on 30 A scales, therefore ignoring 
narrow absorbers does not significantly impact the pho- 
tometry produced here. The instrumental background is 
computed in a similar manner, with extraction boxes off- 
set below the active science region of the detector. The 
instrumental background contributes < 0.6% of the total 
photons measured in each extraction region. 

So that photometry from the blue G130M and red 
G160M gratings can be compared, flux-callibrated spec- 
tra for each exposure (~ 20 — 28 min) are fitted with 
a power law of the form /(A) = f(Xo) X (A/Ao)" A nor- 
malized at Ao = 1400 A, a wavelength covered by both 
gratings. Raw light curves in units of counts s _1 are ex- 
tracted in the wavelengths intervals defined above and 
in 2- min time intervals. Flux calibrated light curves 
at Ao are obtained by dividing the raw light curves by 
exposure-averaged count rates in the same wavelengths 
intervals, and by multiplying them by /(Ao). The re- 
sult is a nearly-continuous, flux-calibrated light curve for 
S5 0716+714 (Figure |2J) at 1400 A, interupted only by 
four periods of Earth occultation. Henceforth, all times 
are quoted in hours starting from MJD 55922 (0:00 UT 
on December 27, 2011). 

3. A FAST FLARE IN THE FAR-ULTRAVIOLET 

Figure [2] shows the strong source variability during 
our five-orbit (7.3 h) HST observing period; the blazar 
brightened by a factor of ~ 40% (0.36 mag) before fad- 
ing to very nearly the initial intensity. Most dramat- 
ically, the source showed a monotonic flux increase of 
0.17 mag during 0.39 h, with the average variability rate 
of —0.45 magh^ 1 during the second orbit. During the 
third orbit, we observed a quasi- monotonic flux decrease 
of +0.11 mag over 0.26 h, with an average variability 
rate of +0.36 mag h _1 . 

Even though there are gaps in our light curve due to 
Earth occultation, the variability appears to be fairly 
simple. Hence, we model the light curve with a con- 
stant background flux level Fq and a single flare template 
which is described by a function 

2 Fi 

Mt) = cxp(4^)+exp(^) (1) 
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Fig. 2. — Far-UV light curve of S5 0716+714, presented as dereddened XFx values at A = 1400A extracted in 120 s time bins. Blue and red 
points show data collected using the G130M and G160M gratings, respectively. The overall light curve was fitted with two models consisting 
of one and two flare templates. The right axis shows magnitude calculated relative to the flux density level of 7 X 10 — 11 ergs" 1 cm -2 . 



TABLE 1 

Flare Model Parameters 3 



Parameter Model A 



Model B 



Fi 
ti 

T r ,l 

T d)1 
F 2 
ti 

T r ,2 



4.78 ±0.03 
1.37 ±0.07 
3.97 ± 0.05 
0.13 ±0.04 
1.85 ±0.12 



4.79 ±0.02 
1.59 ±0.12 
3.99 ±0.05 
0.12 ±0.03 
1.28 ±0.12 
0.25 ± 0.04 
7.34 ± 0.02 
1.40 ±0.53 
0.02 ± 0.02 
3.2 



a Flux density in 10 erg cm 1 s 1 
and times in hours. Peak times are in 
hours since MJD 55922. 

with four parameters: flare normalization Fi, flux rais- 
ing time scale T r .i, flux decay time scale Td.i, and 
epoch ti correspo nding roughly to the flare maximum 
(|Abdo et al.ll2010[) . This Model A is shown in Figure [2 
and its parameters are reported in Table 1. The fit qual- 
ity is mediocre (x 2 = 5.1) and the largest residuals are 
observed during the late stage of the flare. Hence, we 
consider a more detailed Model B, which includes a sec- 
ond flare template. This model provides a substantial im- 
provement in the fit quality (x 2 = 3.2), and the remain- 
ing residuals do not show any long-term structure. The 
main difference in the main flare parameters in Model B, 
compared to Model A, are higher flare amplitude F± and 
shorter flux decay time scale 7d,i. The background flux 
level Fq, flux raising time scale T r i, and the moment of 
peak flux t\ are consistent between models A and B. The 
second flare template f2(t) of Model B has poorly con- 
strained characteristic time scales; in particular, the very 
fast decay time is almost entirely determined by the 
final 120 s bin at the end of the fourth orbit. In Section 
5, we discuss the physical implications of the main flare 
parameters. 

3.1. Spectral Energy Distribution 

Broad spectral coverage allows us to monitor any 
changes in the spectral index of S5 0716+714. Though 
the spectral "throw" of any individual observation is rel- 



atively small (AA/A ~ 23%), they show significantly dif- 
ferent spectral indices during the course of the HST ob- 
servations. We convert our measured ax values (Sec- 
tion 2) to the more conventional a v via the relation - 
ship a u = -(2 + a x ) (|Shull. Stevans fc Danforthll2012l) . 
The SED softens monatonically at the rate of Aa„ ~ 
—0.05 h _1 from a v « —1.0 at the beginning of the 
observation to a v ~ —1.2 near the flare maximum to 
a v w —1.4 at the end of the observation as shown in 
Fig. [3] This trend is consistent with measured flux ra- 
tios from the blue and red ends of each spectral range. 

Plotting spectral index versus intensity (right panel 
of Fig. [3]) gives what could be interpretted as a clock- 
wise hysteresis loop although a v decreases monatonically 
throughout our observation and the loop does not close. 
iKirk. Rieger. fc Mastichiadisl (|1998| ) model electron ac- 
celeration in relativistic shocks and see that the direc- 
tionality of the loop can be related to the propagation of 
a flare through the SED of the blazar. A clockwise loop 
in a^-flux space is typical for blazars observed at ener- 
gies lower than their synchrotron peak. S5 0716+714 is a 
high-peaked BL Lac (HBL) object peaking in the X-rays, 
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Fig. 3. — Power law fits to the individual COS observations 
show a softening of the SED (decreasing a„) throughout the HST 
observation (left panel). Blue circles show fits to the G130M 
(1140 < A < 1450 A) exposures while red diamonds denote G160M 
(1400 < A < 1795 A) exposures. Dotted blue and red lines show 
the fitted index for coadded G130M and G160M exposures. The 
spectral index versus extrapolated flux at 1400 A (right panel) sug- 
gests a clockwise hysteresis loop. 
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so this is the expected behavior as synchrotron heating 
or cooling information propagates from higher to lower 
energies through the SED. However this is in the oppo- 
site sence to the loop seen in the optical observations of 
(ICarini. Walters, fc Hopped l27)Tl . 

The highest-quality spectra can often be obtained by 
scaling the individual exposures and coadding them onto 
a common wavelength scale. However, the significant 
change in the spectral index of S5 0716+714 during the 
observations makes any SED derived from the combined 
observations suspect. Combining just the observations in 
the same far-UV waveband produce slopes of a u « —1.03 
and -1.26 for G130M and G160M channel coadditions, 
respectively (dotted horizontal lines in Figure [3]). Both 
slope values are typical of those of the individual obser- 
vations, but we caution against trying to determine an 
overall SED from the combined far-UV observations. 

4. REDSHIFT OF S5 0716+714 

Determining the source redshift of featureless AGN can 
be difficult, especially at higher redshift where the host 
galaxy cannot be easily identified. Redshift limits can 
be inferred through indirect means (host galaxy non- 
detection, gamma ray emission, etc.). A direct lower 
limit to the source redshift can be determined through in- 
termediate absorption lines in moderate-resolution data 
(usually in the UV). An automated line-finding algo- 
rithm (Danforth et al. 2013, in prep) finds ~ 37 inter- 
vening Lya absorption systems detected to greater than 
4cr significance in the COS spectrum of S5 0716+714. 
Figure 0] shows the observed COS spectral range with 
measured Lya systems along the top edge in both wave- 
length and redshift space. The reddest of these lines is 
at 1497 A, corresponding to z a bs = 0.2315, is confirmed 
by corresponding detections in higher-order Lyman lines 
(Ly/3, Ly7) as well as metal ions (O VI, N V) thus set- 
ting a firm lower limit on the redshift of S5 0716+714. 
The remaining spectral range (1497A < A < 1795A; 
0.232 < 

z Lya < 0.47) is devoid of Lya absorption to a 4cr 
limiting equivalent width of ~ 30 mA (log Vhi ~ 12 .7). 
This technique was first used in lDanforth et al.l (j2010fl to 
constrain the upper redshift limit of 1ES 1553+113 and 
we refine the methodology here. 

We calculate the expected density of intervening Lya 
systems per unit redshift, dAf(z)/dz (shaded curve in 
Figure 2]) . The minimum column density absorber that 
can be detected at a 4<r level is predic ted based on the 
S/N in the data (|Keenev et all I2012D and the dM/dz 
at that limiting column is drawn from the large, low- 
redshift IGM surveys of IDanforth fe Shulll ((20081) and 
ITilton et al.l(l20Tl . We assume no evolution in the Lya 
forest over this small redshift, but including a modest 
evolution produces negligible changes in the predicted 
absorber distribution. The overall density of absorbers 
toward S5 0716+714 (black data points) follows the pre- 
dicted curve consistently though with variations in indi- 
vidual redshift bins most likely due to cosmic variance. 

To constrain the upper limit on source redshift, we 
truncate the predicted dM/dz curve at a range of red- 
shifts and compare the modeled absorber distribution to 
the observed distribution using a K-S test. S5 0716+714 
can be constrained to 0.2315 < z < 0.2668 (68% con- 
fidence), 0.2315 < z < 0.3035 (95%), or 0.2315 < z < 
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Fig. 4. — Intervening absorption lines provide the source redshift 
of S5 0716+714. Vertical ticks at the top of the figure show the 
wavelength/redshift location of all ~ 37 intervening Lyo systems 
seen towards the blazar. The line density as a function of redshift, 
dj\f(z)/dz is shown in bins of dz = 0.03 and compared with the 
expected line density (gray) based on the quality of the COS data 
and the distribution of absorbers at low-redshift in a statistical 
sample HDanforth fc Shullll200a ITilton et al.1120121) . We truncate 
the statistical sample at different redshifts and compare to the ob- 
served distribution with a K-S test. The upper limit is constrained 
to z < 0.304 (z < 0.341) to a 2<r (3a confidence limit). The strong 
absorption system at z = 0.2315 defines the lower redshift limit. 

0.3407 (99%). 

Weak Lya emission is occasionally seen from B LLac 
objects (jStocke. Danforth. fe Perlman et all 1201 lh and 
is an unambiguous, direct way of determining the red- 
shift of blazars to high precision. The only devia- 
tions from a smooth power-law continuum seen in the 
coadded G160M data are likely artifacts produced by 
the coadding individual exposures with rapidly changing 
spectral index and none is convincing as Lya emission 
at z > 0.2315. We set a 3<r upper limit of /(Lya) < 
5 x 10~ 15 erg cm" 2 s _1 (EW < 115 m A) by measur- 
ing the equivalent width detection limit (|Keenev et al.l 
2012) for an emission feature with FWHM« 3 A in 
the rest frame (typical of those seen by Stocke et al. in 
Mrk421, PKS 2005-489 and Mrk501). This intensity 
limit translates to an isotropic luminosity of L(Lya) < 
6.5 x 10 41 erg s _1 . This upper limit is several times 
greater the Lya luminosity observed in lower-redshift 
BL Lacs, so it is not surprising that no emission is seen 
in S5 0716+714. 

5. DISCUSSION 

We performed the deepest and most precise far-UV 
spectroscopic observation of S5 0716+714 to date. We 
did not identify any emission lines, and thus the source 
redshift remains unknown. However, the detection of 
narrow Lya forest features enables the deermination of 
a new, independent constraint — 0.2315 < z < 0.304 at 
the 2a confidence level. This redshift r ange is consistent 
with the result of iNilsson etafl (HI), z = 0.31 ± 0.08, 
based on the photometric detection of the host galaxy 
as well as the intriguing possibility that S5 0716+714 is 
associated with a group of nearby g alaxies observed at 
redshift z = 0.26 IjStickel et al.lll993[i . We note that the 
uncertainty in our estimate is based on the statistical 
properties of the low-redshift Lya forest along many sight 
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lines, while the uncertainty of the iNilsson et al.l (|2008f ) 
estimate is based on the statistics of the host galaxy lu- 
minosities for BLLacs. Our source redshift estimate is 
also consiste nt with the crud e upper limit z < 0.21 ±0.09 
obtained by iPrandini et al.1 <j2010f ) from modeling the 
source spectral energy distribution at V HE gamma-rays 
and a lso with the z < 0.5 estimate of (|Anderhub et al.l 
f2009h . 

Our observations of S5 0716+714 coincided with a 
period of rapid flux variability. This is not surpris- 
ing, since this source is well known for almost uninter- 
rupted variability at multiple time scales. We detected 
episodes of very fast variability rate, up to 0.45 magh -1 . 
This higher than the maximum observed variability rate 
found in the most ext ensive study of IDV in this source 
(jMontagni et al.l 120061). but comparab le to the variabil- 
ity rate found by[Chandra et al. (2011) during one night 
out of five. It seems that such fast variability rates are 
limited to the periods of a fraction of hour. This ap- 
pears to be consistent with the general trend that faster 
variabili ty rates last for short er time intervals (see Fig- 
ure 7 in lMontagni et aLll2006D . It is also known that the 
frequency of light curve segments with constant variabil- 
ity rate decreases exponentially with increasing variabil- 
ity rate. This would indicate rapid flares such as the 
one we observed with HST/COS, while spectacular, are 
probably not importa nt in the overall source energetics. 
However, the study of IMontagni et all (|2006D uses much 
longer sampling rates, and to obtain a statistical picture 
of intra-hour variability of this blazar one needs a much 
larger sample of high-cadence observations. 

The decomposition of the observed light curve of 
S5 0716+714 into a coherent flaring component and a 
quasi-static background (Model B, above) allows us to 
constrain the physics of the emitting region producing 
the flare. The flux raising time scale of the flaring com- 
ponent, T r .i = (7.4 ± 1.8) min, indicates the emitting 
region radius 

R < VcT T>1 /{l + z) ~ 3.4 x 10~ 5 (r>/10) pc (2) 

(for z = 0.26), where T> is the Doppler factor of the emit- 
ting region. The flux decay time constrains the cooling 
time scale 



-^d,l ^ ^cool.obs — 



3(1 + z)m e c 
4£>cr T '"o7 c 



(3) 



where 



l rad 



(1 + q)B' 2 /{&n) is the co- 



moving energy density of the magnetic field (determin- 
ing the synchrotron cooling rate) plus the diffuse radi- 
ation (determining the inverse-Compton cooling rate), 
q ~ L\c/L sya is the Compton dominance parameter 
(which for BL Lacs is of order of unity), and 7 is the 
characteristic random Lorentz factor of electrons produc- 
ing the observed FUV emission. Since the FUV contin- 
uum emission is produced by the synchrotron process, 
we can use the expression for the observed frequency of 
the synchrotron radiation, 



syn,obs 



0.274 Ve^B' ^ c 
(1 + z)m c c ~ 1400 A' 



(4) 



field strength in the emitting region 

B' > 2 G x (P/10)- 1/3 (l + q)- 2/3 . 



(5) 



It follows that 

7 C < 5300 x (£>/10)~ 1/3 (l + q) 1/3 . (6) 
The apparent flare luminosity is 

Li = AirdlFi ~ 3.2 x 10 45 ergs" 1 . (7) 

We use it to estimate the co-moving energy density of 
the synchrotron radiation 



< yn = i!/(4^cP 4 i? 2 ) 

>78 erg cm -3 x (X>/10)~ 



(8) 



In a typical BL Lac, one can neglect the external radi- 
ation, and thus the total energy density of the diffuse 
radiation is u' lad ~ u' syn . Using the estimate for B', we 
can write a direct relation between the Compton domi- 
nance parameter and the Doppler factor: 



q/(l + q) 4/3 ~ 485(D/10)- 16/3 . 



(9) 



We can further use the apparent flare luminosity to es- 
timate the number of electrons contributing to the FUV 
emission 

N^&TiLi/iaTcD^B' 2 ^) (10) 
-2.7 x 10 48 x (P/10)- 8/3 (l + q) 2/3 , 
and their co-moving energy density 

u' c ~ 37V 7 c m c 2 /(47ri? 3 ) 

-2400 erg cm" 3 x (£>/10)~ 6 (l + g). 

If we assume equipartit ion between electron s and the 
magnetic field, u' e 



(11) 



(|Bottcher et all [20091. we will 
obtain a very large Doppler factor, T> ~ 62, and a very 
small Compton dominance parameter q ~ 0.031. Such 
a scenario poses a severe problem of explaining the ex- 
tremely efficient acceleration of the emitting region, but 
it predicts a negligible gamma-ray signature of the FUV 
flare via the Synchrotron Self-Compton (SSC) mecha- 
nism. If we instead assume that q = 1, we will still 
require a large Doppler factor, V ~ 38, but also a strong 
departure from the equipartition, with u' c /u' B ~ 62. It 
should be noted that gamm a-ray emission was de tected 
in S5 0716+714 by MAGIC (|Anderhub et al.l l2009). AG- 
ILE (jVittorini et all 120091) and Fermi at the level of 
Lj — (0.5 — 5) x 10 46 ergs -1 . Hence, the Compton domi- 
nance parameter q cannot be much larger than unity and 
some extreme physical parameter of the emitting region 
producing the fast FUV flare is unavoidable. This makes 
our case of a fast FUV flare similar to the fast gamma- 
ray flares observed in other B LLacs (PKS 2155— 304, 
lAharonian et al.ll2007t Mrk 501. lAlbert et al.ll2007l) and 
one FSRQ (PKS 1222+216, 1 Aleksic et al.l 120 1 ID . Thev 
also seem to require extremely high Doppler factors 
(jBegelman et al.l [2008). and a number of theoretical 



scena r ios were proposed to e x plain them dLevinson 
20071 iStern & Poutanenl l2008t iGhisellini fc Tavecchic 



20081 iGiannios et all l2009t iLvutikov fc Lister! 120101: 
Nalewaiko et all 1201 11) . In the case of S5 0716+714, we 



to eliminate 7 e and obtain an estimate of the magnetic cannot constrain the distance of the emitting region from 



() 
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the central black hole, and thus we have greater freedom 
of theoretical scenarios than in the case of PKS 1222+216 
(jTavecchio et al.ll2011h TNalewaiko et al.ll2012l ). 

We observed strong spectral variability with the sys- 
tematic spectral index change of Aa u = —0.4 over 
the course of ~ 7 h. A diagram of spectral index 
vs. flux density shows probable clockwise hysteresis. 
Strong intraday spectral variations in the optical/UV 
band is rarely reported. For example, IWu et all ()2012f ) 
found the B — R color to vary with the amplitude of 
\A(B — R)\ ~ 0.08 within an hour. This corresponds to 
\Aa v \ = \A(B - R) |/0.413 ~ 0.2, which means a signifi- 
cantly faster variability rate. However, we have found no 
case of a systematic spectral index variation over several 
hours. 

Spectroscopic observations of blazars are proving to 
be versatile tool with which to study a range of astro- 
physical processes from galactic structure to cosmology 
to the study of AGN processes. Our observations of 
S5 0716+714 illustrate the power of time-domain sen- 
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sitivity as well. Not only is the Cosmic Origins Spec- 
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ment ever flown, but the temporal resolution can be ex- 
ploited to look for microvariability on very short time 
scales. A number of other blazars have been observed 
with HST/COS as part of a larger program on low- 
rcdshift cosmology. Many cover only a short time period 
(1-3 HST orbits), but those with similar 5-orbit observa- 
tions do not show appreciable variability. We will con- 
tinue to both look for notable variability in AGN obser- 
vations and to determine the source redshifts of poorly- 
constrained blazars through spectroscopy of their inter- 
vening absorbers. 
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